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We study the time evolution of a system of interacting bosons in a harmonic trap. In the 
low-energy regime, the quantum system is not ergodic and displays rather large fluctuations of the 
ground state occupation number. In the high energy regime of classical physics we find nonergodic 
behavior for modest numbers of trapped particles. We give two conditions that assure the ergodic 
behavior of the quantum system even below the condensation temperature. 
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I. INTRODUCTION 



The understanding of the formation and growth of atomic Bose Einstein condensates is of considerable interest 
and has been investigated experimentally and theoretically . Bose Einstein condensates may be formed while the 
trapped bosons are in contact with a heat bath and a particle reservoir (e.g. via evaporative cooling or from the 
evolution of a nonequilibrium state in an isolated system. In a recent experiment at the MIT JTT1 , harmonically trapped 
bosons were rapidly cooled below the transition temperature, and the subsequent relaxation towards a Bose Einstein 
condensate was observed and measured. The existing theory uses the Boltzmann equation, which is fundamentally 
based on the ergodic StoBzahlansatz. However, this assumption deserves closer scrutiny in the context of bosons in 
, traps. In particular, the motion of two particles confined to a harmonic trap is integrable for any interaction potential 
that depends only on the distance. Closely related to ergodicity is the question about the fluctuations of the ground 
£f} \ state occupation number in an isolated interacting system. For the non-interacting system, these fluctuations have 
t — ■ been calculated recently |l2|,|l3|]. Fluctuations of the interacting system are first computed here. 

The paper is organized as follows. The low-energy quantum Hamiltonian and observables of interest are introduced 
in the second section. In the third section, we present numerical results of the time evolution of the ground state 
occupation number and its fluctuations for bosons in a harmonic and a square well potential, respectively. The 
comparison with a chaotic and ergodic random matrix model is presented in section four. The structure of classical 
phase space for a system of harmonically trapped hard-sphere particles is discussed in the fifth section. Finally, the 
results are summarized in the conclusion. 



C . 

O . II. HAMILTONIAN AND OBSERVABLES 

O 

In this section we describe the quantum Hamiltonian and the observables that may indicate nonergodic behavior. 
We will consider interacting bosons confined to a harmonic potential or a square well potential. Magnetic traps, as 
used in recent experiments, are very well approximated by harmonic potentials. However, the square well potential 
is also of considerable interest, since many theoretical results were obtained for such traps and since the underlying 
classical system of hard-sphere bosons is chaotic and ergodic Q . 

Independent of the specific trap potential, the Hamiltonian may be written 



Here 



is the one-body trap Hamiltonian and 



H = H + V. (1) 

i?o = X>' a K' ( 2 ) 
j 

V = A ^ Vijia ala^auai. (3) 



the two-body interaction. The operators a,j and at annihilate and create one boson in the single-particle trap state 
\j) with energy Ej, respectively. They fulfill the usual bosonic commutation rules. In what follows, the interaction is 
chosen to be a contact interaction. 
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We will study the quantum dynamics of a TV-body system of bosons using the occupation number representation. 
The basis of states is \a) = \no, n\, . . . , nu) with J2i=o k n i ~ N an< ^ "j%l n 0j n i, ■ ■ ■ , n k) — n 3 -\no, n\, . . . , rife). Here nj 
denotes the occupation of the j th single particle state Obviously, the trap Hamiltonian is diagonal in this basis. 

We also note that condensation occurs in d-dimensional harmonic traps for energies E k, N 1+1 / d huj and this 

is the regime we are most interested in. No Bose Einstein condensation occurs in square well potentials in two spatial 
dimensions. 



A. Harmonic Trap 



To be specific, we now consider the case of an isotropic harmonic trap in two spatial dimensions. The single particle 
states are \j) = \rij, rrij) with energy Ej — rijTiuj and angular momentum rrijh {rtij = —nj, —nj + 2, . . . , rij — 2, nj). 
We denote the oscillator wave functions as 4>j(x) = The ground state energy is set to zero. In terms of the 

dimensionless coordinate x — f/do (where do = \fhjmuj sets the length scale of the trap), the contact interaction (||) 
has matrix elements 



V^m = J d x^l (£)(p*j(x)(pk{x)(()i(x) 

s-m fc +ni( 



} ( dXX Af„, + 1 \ mi \ (x) Mnj + l \rn. s \(x) 

M„ h+1 lmkl (x)Mn l+ l lm t |(g). (4) 



Here M /1 ^ K (x) denotes the Whittaker function |l6[ ] and the <5-function ensures conservation of angular momentum. 
The coupling is given by A = Airtuva/do, where a is the s-wave scattering length. Obviously, the total energy E and 
the total angular momentum M are conserved quantities. Furthermore, the motion of the center of mass decouples 
from the single-particle motion yielding a ladder of states spaced by htu for each spatial dimension |l7j . Pitaevskii 
and Rosch have shown, that Hamiltonian ([j]) possesses an additional SO (2, 1) symmetry pq| . This leads to another 
conserved quantity B that can be associated with breathing modes. For fixed quantum numbers M and B the 
spectrum of Hamiltonian (|l|) consists of an infinite ladder with spacing 2huj. The contact interaction (|||]) is valid in 
a dilute system at low energies yjTiui / E 3> a/do, i-e. the scattering length is small compared to the wavelength. 

One may think of the above model as a realization of a i/iree-dimensional oblate trap with u <C oj z and a/ do <C 
a^/mLo z /h <C 1. Under these conditions, the Hamiltonian ([!]) provides a reliable approximation for a dilute gas of 
atoms confined in an oblate magnetic trap at sufficiently low excitation energy. The trap used by the JILA group is 
closest to such a trap. 

We use in the numerical calculation a basis spanned by the eigenstates |a, E, M) of the trap Hamiltonian ([[]) with 
total energy E < -E max and angular momentum M. (a accounts for the different states with quantum numbers E, M. 
Their number is Q(E, M); in what follows we may suppress the dependence on E, M .) Within this basis, the matrix 
of the interaction (|3() is set up as follows. After fixing the total angular momentum M and the maximal energy E max 
we choose one (arbitrary) iV-boson state \a, E, M) with E < E max . We then act with the interaction (||) onto |a) and 
onto all the states created by this procedure until the space with energies E < E max and fixed angular momentum 
M is exhausted. The resulting matrix is sparse. As one check on the numerics we note that the obtained spectra 
display ladders with spacings 2Hlo or Tilo resulting from the breathing mode or the center of mass motion [ p"7| |, 
respectively. Results obtained after the exact diagonalization of Hamiltonian (0) were found to be in good agreement 
with mean field theory even for rather small systems Jl9| . 

At fixed angular momentum, the spectrum of the trap (||) consists of equidistant shells of highly degenerate states 
with energy E = 2KHlu. The spacing between the shells is 2huj (since even (odd) E max requires even (odd) M). We 
label the different shells by their energy quantum number K. Increasing the coupling A shifts the entire spectrum 
towards higher energies and lifts the degeneracies partially. However, the spectrum still exhibits shell structure for 
not too large couplings and number of particles. 



B. Square Well Trap 



In case of a square well potential of width do, the creation operator dj creates a boson in the single particle state 
■ = \nj X , njy) with energy Ej = Eo{n 2 j X + n 2 j y ). Here, Eq — h 2 /2mdo sets the energy scale. nj X and nj X are the 
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quantum numbers in the x and y-direction, and the wave functions are </>j(x,y) = 2sm(jij X irx) sm(nj y iry), where 
distances < x, y < 1 are measured in units of do- States may be characterized by their behavior under reflection 
at the axes x = 1/2, y = 1/2. (We do not specify the behavior under reflection x <-> y.) The coupling is given 
by A = AnEoa/do in terms of the s-wave scattering length a. In what follows we consider those states which have 
the same symmetry as the ground state. The Hamiltonian matrix is built up as for the harmonic trap. At zero 
coupling A, the spectrum is highly degenerate and thus exhibits shell structure. With increasing coupling, most of 
the degeneracies are lifted, the entire spectrum is shifted towards higher energies, and the shells overlap very soon. 
This behavior is different from the harmonic trap where the shell structure persists even for relatively large values of 
the coupling. 

C. Observables 

Nonergodic behavior manifests itself in the difference between the time average of a dynamical observable and its 
ensemble average. The most interesting observable is the ground state occupation number, no = a ao- This observable 
has been measured, and its microcanonical fluctuations have been computed recently for the non-interacting system 
jj^fl . The time evolution for its expectation value with the trap eigenstate state \a) is 

n^\t) = (a\U\t)fl U(t)\a). (5) 

Here 

U(t) = exp \-iHt/h] (6) 

is the time evolution operator. In practice, we evaluate eq. (||) by applying U(At) with a suitable chosen At onto 
the initial state \a) by using its Taylor expansion []. This procedure allows one to treat rather large matrices. The 
unitarity of U{t) was numerically ensured by an error of the norm {aU* (t)\U {t)a) that was smaller than one percent. 
The results are compared to the equilibrium value 

^ si 

nfi= -^2{a\h \a). (7) 

a=l 

In the case of the harmonic trap potential, this equilibrium value is the ensemble average over all ft(E,M) states 
within one energy shell. Note that the states within this ensemble have different values of the conserved quantity 
B associated with the breathing mode. In the case of the square well potential the ensemble comprises all states 
within an energy layer SE « 5n 2 h 2 /2 vnd\ around E that have the same symmetry as the ground state. SE is 
reasonably smaller than the ground state energy 2Nit 2 Ti 2 /2m<ig even for modest numbers of particles. Note that 
these ensembles differ from the microcanonical ensembles which are determined by energy only and do not depend 
on additional quantum numbers. However, nonergodic behavior in a sector defined by additional quantum numbers 
implies nonergodic behavior in the entire energy shell. 

A second and more quantitative indicator for nonergodic behavior is given by the fluctuations of the ground state 
occupation number. These arc defined by 



^o(A) = - £<*|n |*> 2 - ( ^ £<*M*) ) ' ( g ) 

where \^(E, M)) are the eigenstates of the interacting system. Note that the fluctuations (||) depend on A. 

For an ergodic system, these fluctuations should be very small since an expectation value should not depend on the 
particularly chosen state within an energy shell of Q, states. It is useful to compare the fluctuations for the interacting 
system with the fluctuations for the non-interacting system 

5n 2 (0)=^ 2 -m 2 = ^J2^ a \ fl ^)~mf. (9) 

a 

To compute the fluctuations (||), we diagonalize the Hamiltonian (pi) in the basis \a) and obtain the eigenstates 



1 We have used a step size XcoAt ~ 0(1/10) to assure a fast convergence of the Taylor expansion. 
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III. NUMERICAL RESULTS 



We describe the numerical results obtained for the harmonic trap and the square well potential. 

A. Harmonic Trap 

We choose A = 0.15hu>, N = 10 and E max = 12?iuj, This space has a dimension of 1530 and according to the 
criterion given in Refs. should exhibit condensation. Fig. |l| shows n a \t) for a few states of the (K — 4)-shell 
that have initially 2 < Uq (t — 0) < A — 1. These states have occupation numbers that are initially different from 

the equilibrium value no defined in eq. (0). As shown in Fig. 0, for n ( o a) {0) < n [n "' (0) > n ] the ground state 
occupation number rises [falls] on a time scale involving many periods and saturates on a higher [lower] level. However, 
the saturation is lower [higher] than the equilibrium occupation number and is obviously different for different initial 
states. This indicates that the quantum system is not ergodic. 

Table. | shows the fluctuations (||) for different values of the coupling A and TV = 10 bosons. Again, we have chosen 
E max = 12hui. With increasing A the spectrum is shifted towards higher energies but the shell structure persists. For 
nonzero A the fluctuations are suppressed in comparison to the non-interacting case, and are practically independent 
of A for < A < 0.2. The last row of Table [| presents the fluctuations one would expect in a ergodic system. These 
findings confirm the results obtained from the time evolution. They show that the quantum system is not ergodic in 
the regime of low energies and modest number of bosons. 

We a lso consider the A-dependence of the fluctuations. Fixing A = 0.025^ and £ mal = 12huj, Table. |l| shows the 
ratio y/ Stiq(A) / 5riQ(0) for different numbers of particles. One observes that the fluctuations increase with increasing 
N and approach the fluctuations of the non-interacting system. This can be understood as follows. Since E max has 
been fixed, the energy available per boson is decreasing with increasing numbers of bosons and thus the regime N ^> 1 
corresponds to the low temperature regime where most bosons are in the condensate. The interactions induce mostly 
scattering within the condensate, and this corresponds to the diagonal part of the interaction matrix ([|), which does 
not yield any change in the cigenfunctions. 

As shown in Section [v] the classical system becomes chaotic for A 5 / 2 (a/ do) ^JTioj/E > 1. Thus, the classical system 
is chaotic for the values N = 40, A = 0.025?iaj and E = 6Hlo used in Table [n]. Nevertheless, this property of the 
classical system is not reflected by the quantum system. 

B. Square Well Potential 

We next consider a square well potential, keeping the same number of bosons (A = 10) as before. Fig. || shows 
the fluctuations of the ground state occupation number as a function of energy for different values of the coupling 
A. With increasing values of A the results have been shifted to lower energies, such that the average ground state 
occupation number gets (almost) independent of A. Fig. |^ shows that the fluctuations are relatively large even for 
the interacting system. Thus the quantum system is not ergodic at low energies. We recall that a classical system of 
hard-sphere particles confined to a square well is chaotic and ergodic jl4| . This classical behavior is expected to be 
followed by the quantum system in the semiclassical limit pp| , |2l|] when the wave length is sufficiently small compared 
to other length scales including the scattering length. However, our contact interaction demands the opposite limit 
of long wave lengths. One also recognizes in Fig. [2| that the fluctuations decrease with increasing A. One may only 
speculate whether this behavior is tied to the ergodicity of the classical system. 

Note that nonergodic behavior of a classically chaotic system is not unexpected in the low-energy quantum regime. 
The same observation has been made, e.g. for a system of one particle confined to the stadium billiard [ p2| . 

IV. RANDOM INTERACTIONS 

We would now like to contrast the results presented in the last section with a quantum Hamiltonian that does 
display ergodic behavior. To this purpose, we restrict ourselves to the states within one energy shell and model 
the interaction by a f2(-E, M)-dimensional random matrix which is drawn from the Gaussian Orthogonal Ensemble 
(GOE). Within one shell, the results depend only trivially on A and the variance of the random matrix elements. Our 
random matrix model is motivated by the observation that, within the semiclassical regime, fluctuation properties 
concerning eigenvalues and wave functions of classically chaotic systems are universal and coincide with those of the 
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Gaussian random matrix ensembles |2jjj21|]. Thus, one would expect that a random interaction would yield ergodic 
behavior. The appropriate time evolution of the ground state occupation number is shown in Fig. |^. As expected, 
the equilibrium is practically reached by every initial state, and the fluctuations are relatively small. 

Using random matrix theory, we may also obtain analytical results for fluctuation of the ground state occupation 
number. Let \^f) = J2 a c a\ a ) be an eigenvector of the random matrix Hamiltonian. Within the GOE the coefficients 
c a are uniformly distributed over a Sl-dimensional sphere of unit radius. The normalized joint probability distribution 
for k coefficients is given by pij 



^■■■■■^'^rw^mi 1 -^/") ■ (10) 

To compute expectation values, the average over the energy shell is now replaced by the GOE average, and we denote 
the latter by brackets (.). This is justified for £1 3> 1 since the GOE becomes ergodic in the limit of infinitely many 
levels, i.e. each of its members displays the same fluctuation properties as the ensemble The ground state 

occupation number of the eigenstate |\?) is 



nW = (*K|*) = £c^\ (11) 

a 

and its GOE average is 

l 

= / dcP^n^ =nE (12) 



as expected. Using P2 we obtain for the variance 

5njj(GOE) = ((n™ - (n«)) 2 ) = ^^(0), (13) 

where <5tiq(0) is the variance of the non-interacting system given in eq. @. This shows that random interactions 
lead to a tremendous suppression in comparison to the non interacting case and vanish in the limit of infinitely 



many levels. Note that an expression similar to eq.(13) holds for any one-body observable that commutes with the 
unperturbed Hamiltonian. Table III compares the fluctuations of the ground state occupation number of the non 
interacting system with the numerical results obtained with a random Hamiltonian and with the analytical result jl^). 
Obviously, the fluctuations of the non-interacting system are much larger than those of the random Hamiltonian. 
One also recognizes that the agreement between the analytical result and the numerical simulation improves with 
increasing dimension of the random matrix. These results show that the interaction mediated by s-wave scattering 
is very different from a chaos simulating random many-body interaction. 

V. CLASSICAL PHASE SPACE 

It would be interesting to investigate the ergodic properties also at higher energies. Unfortunately, a treatment of 
the full many-body system is very difficult outside the low-energy quantum regime. However, it is useful to study the 
classical many-body system in more detail. Within the semiclassical regime, where the wave length is small compared 
to any length scale including the scattering length (i.e. (a/dp) ^/ E /hm 1 [|), the quantum system is expected to 
reflect the properties of the underlying classical system |2(],|2lj . Sinai has shown that the classical system of hard- 
sphere particles in a square well potential is ergodic and chaotic jl4| . We therefore consider the classical dynamics of a 
dilute system of hard-sphere particles of radius a confined to a two-dimensional isotropic harmonic oscillator. At low 
energies, the corresponding quantum system is governed by the Hamiltonian ([!]). Of course, in the low-energy limit, 
any short ranged two-body interaction with s-wave scattering length a would yield the quantum Hamiltonian (jl]), 
and there is no reason to choose a hard-sphere interaction. However, hard potentials induce more chaos than softer 



2 This estimates is based on the assumption that one particle has the entire energy E. A more conservative estimate would 
use the fraction E/N for single-particle energies. 
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ones, and they are easier to treat. This justifies and motivates our specific choice. Nevertheless, most of the results 
derived below are valid for any two-body interaction that depends only on the distance and vanishes for distances 
larger than the scattering length. 

Let us first consider the two-body system with coordinates and momenta f\ and pi (i — 1,2), respectively. Intro- 
ducing coordinates and momenta q± = (fi ± r^/v^, p± = (p 1 ±p 2 )/y/2 shows that the the system is equivalent 
to a system of two nonidentical and non-interacting particles in a harmonic oscillator where the particle with coor- 
dinates (f_ 'sees' a spherical symmetric hard core potential with radius \f 7 ia in addition to the confining harmonic 
potential. In the new coordinate system, single particle energies and single particle angular momenta are conserved, 
which renders the system integrable. Obviously, scattering occurs only if the angular momentum of the particle with 
coordinates <f_ is sufficiently low, i.e. < 4h 2 (a/do) 2 [E/huj — (a/do) 2 ]- Let us assume that the initial conditions of 
the two-body system are uniformly distributed over the shell of energy E. The probability of scattering is basically 
given by integrating Q(ATi 2 {a / do) 2 [E / Tiuj — (a/d ) 2 ] — l 2 _) over the energy shell in phase space. One obtains 

m f 1 : 1 < t < 2 

P {2) (0 = < 24/5 9 , 8 . 2<C (14) 



where £ = (do / a) 2 E / {Tllo) is a dimensionless parameter. Note that £ ^> 1 in the traps used in recent experiments. 
Obviously, the scattering probability is very small for two-body systems, and this is a consequence of the well known 
fact that the periods do not depend on the energy in harmonic potentials. For the A-body system we may use eq. 
( |l4| ) to compute the probability of having at least one scattering. We give an estimate. On average, any two particles 
have a fraction 2/N of the total energy. There are N(N — l)/2 different pairs of particles which may collide. Thus, 



p (N )ftj N(N- 1) p(2)(2W 



(15) 



A more accurate calculation of the leading term oc jV 5 / 2 ^ 1 / 2 in eq. ( |l5"l) confirms that the estimate is a good 
approximation. Note that formulae eq. dig ) and eq. (|lj) are valid for any short ranged two-body interaction that 
vanishes for distances larger than a. The absence of any scattering for some fraction of the energy shell is a consequence 
of the energy independence of the periods in the isotropic harmonic potential. Thus, the interacting classical system 
is not ergodic if the number of particles is not too large. 

Let us also consider those regions of phase space that involve scattering among more than two particles. Starting 
trajectories in such regions, we computed the Lyapunov exponent [§5j for systems with 2 < N < 25 particles. As 
expected we found positive Lyapunov exponents. Thus, scattering between more than two particles yields chaotic 
dynamics in the corresponding fraction of phase space. As a rule of thumb, eq. (|l5|) shows that the classical system 



becomes chaotic for N 5 / 2 (a/d )^huj/E > 1. 

Let us also discuss the case of a three-dimensional cylindrical symmetric trap with axial symmetry. In the absence 
of interactions the motion in the z-direction decouples from the motion in the radial plane. Since a scattering in three 
dimensions also is a scattering of the particles in the radial plane, the results ( |l4| , |l5|) derived in this section also hold 
for three-dimensional axially symmetric traps provided the total energy is replaced by the radial energy. 

Our results are consistent with the recent theoretical observation that the quasiparticle motion of collective and 
single-particle excitations of Bose Einstein condensates is only weakly chaotic pq] . 



VI. CONCLUSION 



We have studied the ergodic properties of trapped bosons that interact via s-wave scattering. In the low-energy 
quantum regime we find nonergodic behavior for both the harmonic trap and the square well potential. The nonergodic 
behavior may be seen in the difference between time averages and ensemble averages and in rather large fluctuations 
of the ground state occupation number. For modest numbers of harmonically trapped bosons the fluctuations are 
smaller than for the non-interacting system and almost independent of exact magnitude of the s-wave scattering 
length. For larger numbers of harmonically trapped bosons the fluctuations of the interacting system approach the 
fluctuations of the non-interacting system in the low-temperature limit. This shows that even the many-body system 
is not ergodic at sufficiently low temperatures. 

The analysis of the classical phase space structure shows that, unlike square well potentials, harmonic potentials 
do not necessarily lead to chaotic behavior of trapped interacting particles. Only for N 5 ^ 2 (a/do)^huj/ E > 1 a 
considerable fraction of classical phase space is chaotic. This is achieved with presently used traps if A^ > 40 particles 
are trapped at temperatures of the order of the condensation temperature [E ~ N 1+1 ^ d hcj). The quantum system 
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is expected to follow the classical system in th e semi classical regime where the wave length is small compared to the 
scattering length. This is the case for (a / do) \J E / Tito 3> 1, and N 3 • 10 4 or N 3> 10 4 atoms have to be trapped in 
d = 3 or d = 2 dimensions, respectively. 

The parameter that compares the strength of the interaction with the kinetic energy is Na/do Ml. For the 
experiments with rubidium at JILA [Q, lithium at Rice University and sodium at the MIT |3j, one has 
Na/do < 1, Na/do > 1, and Na/do 3> 1, respectively. Our computations are restricted to the weakly interact- 
ing regime Na/do < 1, and as one consequence we do observe shell structure. It is not obvi ous, that this restriction 
is also responsible for the nonergodic behavior observed in this work. The results of Section IV show that any small 
random interaction may render the system ergodic since the noninteracting system is highly degenerate. Whether the 
quantum system is ergodic also in the low energy quantum regime for Na/do > 1 remains an open question. Given the 
fluctuations ( p"3| ) for the system with random interactions on the one hand and the results p"2| ] for the non-interacting 
system on the other hand, experiments should be able to reveal how ergodic systems of trapped atoms really are. 

We have considered the fluctuation of the ground state occupation number as one observable that is sensitive 
to nonergodic behavior. Chaotic systems exhibit level repulsion within sectors of definite symmetry and may be 
distinguished from integrable ones by their level statistics [P3| , ^l| . Due to its 5*0(2, 1) symmetry the two-dimensional 
harmonic trap with contact interaction is quite special. Sectors of definite symmetry (fixed angular momentum M 
and fixed value of B) are ladders of levels with spacing 2huj |H and are expected to display the same level statistics 
as a harmonic oscillator. 

It is also interesting to discuss the growth of the condensate in harmonic traps. As long as the shells do not overlap 
and first order perturbation theory in A is valid, the quantum mechanical time scale set by the interaction is given 
by the product 1/Xlo. This may be compared to an approach using the Boltzmann equation, where the time scale is 
proportional to 1/A 2 (when transition rates are obtained from Fermi's golden rule). Thus, for sufficiently small values 
of A, the use of the Boltzmann equation in combination with transition rates resulting from Fermi's golden rule is 
inappropriate. It yields times for the condensate formation that are too large. This finding is a consequence of the 
high degeneracy of the harmonic trap spectrum. 



ACKNOWLEDGEMENT 



We acknowledge discussions with Paulo Bedaque and Aurel Bulgac. This work was supported by the Department 
of Energy on Contract No. DE-FG-06-90ER-40561. 



M.N. Anderson, J.R. Ensher, M.R. Matthews, C.E. Wieman, E.A. Cornell, Science 269 (1995) 198; 
C.C. Bradley, C.A. Sacket, J.J. Toilet, R.G. Hulet, Phys. Rev. Lett. 75 (1995) 1687; 

K.B. Davis, M.-O. Mewes, M.R. Andrews, N.J. van Druten, D.S. Durfee, D.M. Kurn, W. Ketterle, Phys. Rev. Lett. 75 
(1995) 3969 

cond-ma t/9805034 to be published in Rev. Mod. Phys. 



For a review see, e.g. F. Dalfovo et al, 



H.T.C. Stoof, :ond-mat/9805393, to be published in J. Low Temp. Phys. 



M. Holland, J. Williams, J. Cooper, Phys. Rev. A 55 (1997) 3670; 

H. Wu, E. Arimondo, C.J. Foot, Phys. Rev. A 56 (1997) 560; 

D. Jaksch, C.W. Gardiner, P. Zoller, Phys. Rev. A 56 (1997) 575; 



C.W. Gardiner et al, Phys. Rev. Lett. 79 (1997) 1793; |gond-mat/9801027 ; 
For a review see, e.g. W. Ketterle, N.J. van Druten, Adv. Atomic, Mol. and Opt. Phys. 37 (1996) 181 
H.-J. Miesner et al, Science 279 (1998) 1005 
S. Grossmann, M. Holthaus, Phys. Rev. Lett. 79 (1997) 3557; 
M. Gajda, K. Rzazewski, Phys. Rev. Lett. 78 (1997) 2686; 
Ya.G Sinai, Sov. Math. Dokl. 4 (1963) 1818 

SR. de Groot, G.J. Hooman, C.A. Ten Seldam, Proc. Roy. Soc. (London) A203 (1950) 266; 
V. Bagnato, D.E. Pritchard, D. Kleppner, Phys. Rev. A 35 (1987) 4354 

M. Abramowitz, LA. Stegun, Handbook of Mathematical Functions, Dover Pub., New York, 1972 
AL. Fetter, D. Rokhsar, Phys. Rev. A 57 (1998) 1191 
LP. Pitaevskii, A. Rosch, Phys. Rev. A 55 (1998) R853 
T. Haugset, H. Haugerud, Phys. Rev. A 57 3809 

M.-J. Giannoni it et al (Eds.), Chaos in Quantum Physics (North-Holland, Amsterdam, 1991) 



7 



[21] T. Guhr, A. Miiller-Groeling, H.A. Weidenmiiller, Phys. Rep. 299 (1998) 189 

[22] F. Borgonovi. G. Casati, B. Li, Phys. Rev. Lett. 77 (1996) 4744 

[23] O. Bohigas, M.-J. Giannoni, C. Schmit, Phys. Rev. Lett 52 (1984) 1 

[24] T.A. Brody et al, Rev. Mod. Phys. 53 (1981) 385 

[25] G. Benettin, L. Galgani, J.M. Strelcyn, Phys. Rev. A 14 (1976) 2338 

[26] M. Fliesser, A. Csordas, R. Graham, Phys. Rev. A 56 (1997) 4879; 



M. Fliesser, R. Graham, cond-mat/9806115 





K = 2 


K = 3 


K = 4 


K = 5 


A = Q.Ohuj 


0.83 


1.11 


1.32 


1.51 


A = 0.025ftaj 


0.56 


0.55 


0.59 


0.64 


A = 0.05fi,w 


0.54 


0.55 


0.60 


0.64 


A = 0.075?kj 


0.54 


0.55 


0.60 


0.65 


A = Q.Yhw 


0.53 


0.55 


0.59 


0.64 


A = 0.15fi,w 


0.53 


0.55 


0.59 


0.63 


A = Q.2huj 


0.52 


0.55 


0.58 


0.64 


GOE 


0.35 


0.27 


0.18 


0.12 



TABLE I. Fluctuations of the ground state occupation number */ Srig(X) for different shells (labeled by K) and different 
values of the coupling A. Systems of N = 10 bosons are considered. For A = 0, the states of the Kth shell have the energy 
E = 2Khuj. The shells comprise Q(E = 2KHuj,M = 0) = 9,31, 109,339 states for K = 2,3,4,5, respectively. The last row 
shows results expected for an ergodic system. 





K = 2 


K = 3 


K = 4 


K = 5 


N = 12 


0.71 


0.56 


0.51 


0.46 


N = 15 


0.77 


0.65 


0.58 


0.54 


N = 20 


0.83 


0.73 


0.68 


0.63 


N = 25 


0.86 


0.78 


0.73 


0.69 


N = 40 


0.91 


0.85 







TABLE II. Fluctuations of the ground state occupation number y/ Stiq(X) normalized by the fluctuations of the 

non-interacting system yj <5?1q(0) for different shells (labeled by K) and different numbers N of bosons. The coupling is 
fixed to A = 0.025fto;. The shells comprise Q(E = 2K%u, M = 0) = 9, 31, 109, 339 states for K = 2, 3, 4, 5, respectively. 





K = A 


K = 5 


K = 6 


HO 


1.32 


1.51 


1.64 


RM 


0.20 


0.112 


0.072 


GOE 


0.18 


0.116 


0.072 



TABLE III. Fluctuations of the ground state occupation number in different shells labeled by K. Results are shown for 



the harmonic oscillator (HO), harmonic oscillator with random interaction (RM), and the analytical result ( p_3[ ) derived from 
random matrix theory (GOE). The shells comprise f2 = 109, 339, 1039 states for K — 4, 5, 6, respectively. 
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FIG. 1. Time evolution (thin lines) of the ground state occupation number (g) for different initial states of the (K = 4)-shell 
that are initially away from the equilibrium value (thick line). Results are obtained for a system of N = 10 harmonically 
trapped bosons that interact via s-wave scattering. The coupling is A = 0.15hui. 
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FIG. 2. Fluctuations of the ground state occupation number y/ 5n 2 (X) for a system of N = 10 bosons confined to a square 
well and different values of the coupling A. (X/E = 0, ^, ^j, ^, ^, ^ from top to bottom.) The energy is given in 

units of Eo = Tv 2 h 2 /2md,Q. 
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FIG. 3. Time evolution (thin lines) of the ground state occupation number (g) for different initial states of the (K = 4)-shell 
that are initially away from the equilibrium value (thick line). Results are obtained for a system of N = 10 harmonically 
trapped bosons that interact via random interactions. 
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